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ABSTRACT 


Accurate  predictions  and  understanding  of  littoral  and 
coastal  wave  conditions  are  of  major  importance  to  military 
operations  and  civilian  coastal  zone  management.  Although 
WaveWatchlll  (WW3)  is  used  by  many  operational  forecasting 
centers  around  the  world,  there  is  a  lack  of  field  studies 
to  evaluate  its  accuracy  in  regional  applications  and  under 
extreme  conditions,  such  as  hurricanes. 

Data  from  seven  National  Data  Buoy  Center  (NDBC)  buoys 
in  the  Gulf  of  Mexico,  together  with  an  array  of  pressure 
and  pressure-velocity  sensors  deployed  on  the  Florida 
Panhandle  shelf  during  the  Office  of  Naval  Research  (ONR) 
SAX04  experiment,  were  used  to  test  WW3  predictions  of 
extreme  waves  generated  by  Hurricane  Ivan. 

The  model  predicts  large  differences  between  wave 
conditions  on  the  left  and  right  sides  of  the  hurricane 
track  owing  to  the  difference  in  "dwell  time"  between  wave 
propagating  against  and  with  the  storm.  Analysis  reveals  a 
tendency  to  predict  smaller  wave  heights  and  later  arrival 
of  hurricane  swell  than  is  observed.  Additionally,  the 
default  operational  setting  for  dissipation  by  bottom 
friction  yields  too  much  dissipation  on  the  continental 
shelf.  Overall,  the  agreement  between  observations  and 
model  predictions  is  reasonable. 
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I. 


INTRODUCTION 


Since  ancient  times,  when  men  decided  to  explore  the 
unknown  sea,  wind  and  waves  were  always  feared  and  were 
objects  of  major  interest.  In  1490,  the  great  navigator. 
Count  Vasco  da  Gama,  used  to  spend  a  significant  amount  of 
time  per  day  watching  the  sea  with  the  purpose  of 
understanding  its  mysteries,  how  waves  react  to  wind  or  the 
tide  effects.  Five  hundred  years  later  we  still  spend  a 
considerable  amount  of  time  watching  and  monitoring  the  sea 
in  quest  of  the  same  knowledge. 

Due  to  the  present  demands  of  our  society,  we  can  not 
just  understand  its  mechanisms.  We  also  have  to  be  able  to 
predict  it.  In  order  to  carry  out  that  task,  in  the  last 
decades  sophisticated  numerical  models  have  been  developed 
and  implemented  on  high  performance  computers  to  provide 
operational  global  and  regional  wave  forecasts. 

A.  BACKGROUND  AND  MOTIVATION 

Forecast  models  had  their  beginning  with  Sverdrup  and 
Munk  during  the  1940'' s,  based  on  a  semi-empirical  wave 
forecasting  relationship  that  made  extensive  use  of  graphic 
solutions  (Sverdrup  and  Munk,  1947),  but  it  was  Gelci  in 
the  late  1950's  who  developed  the  first  wave  prediction 
model  using  a  balance  equation  between  wind  input  and 
dissipation.  The  second  generation  wave  models  started  in 
the  early  1970' s  with  the  Joint  North  Sea  Wave  Project 
(JONSWAP)  (Hasselmann  et  al . ,  1973) .  Although  an 
improvement  in  forecasting  capabilities  was  achieved,  this 
second  generation  model  still  had  a  simplified  view  of 
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nonlinear  transfer,  and  did  not  handle  the  transition 
between  windsea  and  swell  properly.  In  1985,  the  SWAMP 
study  concluded  that  none  of  the  existing  models  were 
reliable  in  extreme  situations,  and  in  fact,  it  was 
proposed  that  a  third  generation  model  should  be  developed. 

Improvements  in  computational  capabilities  and 
numerical  methods  led  to  the  first  third  generation  wave 
model  WAM  (WAMDI,  1988)  .  WAM  solves  the  full  spectral 
energy  balance  including  nonlinear  wave-wave  interactions. 
The  source  function  was  still  represented  as  a  combination 
of  wind  input,  nonlinear  transfer  and  dissipation,  however 
significant  improvements  were  introduced  in  the 
parameterization  of  the  source  components. 

WAVEWATCHI I I  (WW3) (Tolman,  2003),  another  third 
generation  wave  model,  introduced  further  refinements 
product  based  on  the  previous  models  WaveWatch  I  and  II, 
including  a  third-order  numerical  scheme  and  new 
formulations  of  the  wind  input  and  dissipation.  WW3  is  the 
operational  model  used  at  FNMOC  for  global  and  regional 
wave  forecasts. 

B .  STUDY  PURPOSE 

Natural  disasters  associated  with  severe 
meteorological  disturbances  are  among  the  most  impressive 
and  destructives  events.  In  addition  to  damage  from 
hurricane-force  winds,  extreme  waves  cause  beach  erosion 
and  property  damage.  Accurate  prediction  of  waves  and 
associated  sediment  transport  is  critical  for  coastal  zone 
management . 
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also 


affect 


Large  waves  also  affect  military  operations 
includeding  amphibious  landings,  mine  warfare  and  military 
insertions.  The  extent,  to  which  waves  influence  so  many 
processes  and  operations  in  the  marine  environment,  has 
driven  a  need  for  accurate  wave  prediction  models. 
Presently,  the  main  forecasting  centers  are  using  third 
generation  global  wave  models,  such  as  WW3 .  Numerous 
studies  have  been  performed  to  validate  model  predictions 
with  in  situ  and  remote  sensing  observations,  but  the 
accuracy  of  these  models  in  extreme  hurricane  conditions  is 
not  well  understood. 

The  operational  wave  prediction  model  WW3  was  designed 
for  global  applications  with  relatively  coarse  resolution 
(typically  around  1  degree) .  Its  accuracy  and  stability  has 
been  demonstrated  through  extensive  validation  with  buoy 
and  satellite  altimeter  measurements  (Tolman,  2003) .  The 
accuracy  of  WW3  in  regional  applications,  that  requires 
higher  resolution  to  resolve  coastline  and  bathymetry 


effects,  and  the  performance  of  WW3 
meteorological  conditions,  like  hurricanes, 
known  and  is  the  topic  of  this  study. 


in  extreme 
is  not  well 
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II.  DATA  SETS  AND  METHODOLOGY 


This  study  focuses  on  a  period  of  4  days,  September  13 
to  16,  2004,  when  the  category  V  hurricane  Ivan  crossed  the 
Gulf  of  Mexico  and  made  landfall  on  the  Alabama  coast. 
Fortunately,  at  that  time  the  Sediment  Acoustics  experiment 
funded  by  the  Office  of  Naval  Research  (ONR)  took  place  on 
the  Panhandle  coast,  providing  unique  observations  of  the 
transformation  of  extreme  hurricane-generated  waves  across 
the  continental  shelf.  These  data  were  analyzed  together 
with  available  data  from  permanent  buoys  in  the  Gulf  of 
Mexico  and  the  Geosat  Follow  On  satellite.  The  data  sets 
and  analysis  methodology  are  described  in  this  chapter. 

A.  BUOY  DATA 


100  00' W  95  00' W  90  00' W  85  00'  W  80  00' W 


Figure  1.  The  Hurricane  Ivan  track  (blue  dashed  line)  and 
locations  of  the  seven  NDBC  buoys  that  measured  wind  and 

wave  conditions. 
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The  buoy  datasets  used  in  this  thesis  were  provided  by 
the  National  Data  Buoy  Center  (NDBC)  .  Wind  and  wave  data 
were  collected  by  seven  buoys  in  the  Gulf  of  Mexico.  The 
main  characteristics  and  geographical  positions  of  these 
buoys  are  listed  in  Table  1.  The  buoys  are  deployed  in 
depths  ranging  from  14  to  3200  meters  (Table  1)  .  Buoys 
42003  and  42039  were  very  close  to  the  hurricane  path  and 
recorded  extreme  wind  and  wave  conditions.  Buoy  42040, 
located  along  the  path  was  damaged  during  the  passage  of 
the  hurricane  and  stopped  collecting  data. 


Table  1 .  NDBC  buoys  used  in  the  study 


Buoy 

Position  (WGS-84) 

Lat  (N) /Lon  (W) 

Dir 

Type 

Diameter 

(m) 

Depth 

(m) 

42001 

25.90 

89.67 

Yes 

Discus 

12 

3246 

42003 

26.01 

85.91 

Yes 

Nomad 

6 

3233 

42007 

30.09 

88.77 

Yes 

Discus 

3 

14 

42036 

28.50 

84.52 

Yes 

Discus 

3 

54.5 

42039 

28.79 

86.02 

Yes 

Discus 

3 

291.4 

42040 

29.18 

88.21 

Yes 

Discus 

3 

443.6 

42041 

27.50 

90.46 

Yes 

Discus 

3 

1435 

All  NDBC  buoys  are  equipped  with  a  suite  of 
commercially  available  sensors  such  as  anemometers  and 
barometers  that  provide  hourly  standard  meteorological 
observations  including  mean  wind  speed  and  direction, 
maximum  gust,  barometric  pressure,  air  temperature,  and  sea 
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surface  temperature.  To  collect  wave  data,  the  buoys  are 
equipped  with  accelerometers  that  are  used  to  measure 
vertical  acceleration  with  a  sample  frequency  of  1 . 5Hz  for 
20  minutes.  The  acceleration  time  series  are  inteqrated 
twice  to  form  sea  surface  elevation  time  series 
(acceleration  back  to  velocity  back  to  displacement) .  A 
Fast  Fourier  Transform  (FFT)  is  applied  to  the  data  by  the 
processor  on-board  the  buoy  to  transform  the  data  from  the 
temporal  domain  into  the  frequency  domain  (Steele,  1993)  . 
Spectral  analysis  yields  wave  frequency  spectra  and 
standard  parameters  such  as  siqnificant  wave  heiqht  (SWH) 
and  peak  period  (Tp) .  The  spectral  data  (enerqy  density  and 
directional  moments)  are  computed  for  47  frequency  bands, 
in  the  range  of  0.020Hz  to  0.485Hz  (Table  2) . 


Table  2 .  Frequencies  and  bandwidths  used  in  spectral 

analysis  for  NDBC  buoys 


Frequency  (Hz) 

Bandwidth  (Hz) 

0.020  -  0.100 

0.005 

0.100  -  0.350 

0.010 

0.350  -  0.485 

0.015 

B.  SEDIMENT  ACOUSTICS  EXPERIMENT  2004  (SAX04)  -  IN  SITU 
INSTRUMENTATION 

SAX04,  an  ONR  experiment  aimed  at  understanding 
acoustic  propagation  through  sediments,  was  conducted  off 
Panama  City,  Florida  from  early  September  through  the 
middle  of  November  of  2004.  In  addition  to  an  extensive 
suite  of  acoustic/geological  measurements,  13  wave¬ 
measuring  instruments  were  deployed  by  the  Naval 
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Postgraduate  School  at  9  different  sites  (Figure  2),  and 
these  data  are  used  in  this  study.  The  instrumentation 
includes  4  Nortek  Vector  puv  sensors  (sensors  that  measure 
pressure  and  horizontal  velocity,  u  and  v)  (Figure  3)  ,  8 

pressure  sensors  (Figure  4),  and  1  wavebuoy  (Datawell 
Directional  Waverider) . 


Table  3.  Geographical  position,  depth  and  type  of  sensors 

deployed  at  9  different  sites 


Site 

Latitude  (N) 
(WGS-84) 

Longitude (W) 
(WGS-84) 

Depth (m) 

Sensor 

1 

30°  23.395' 

86°  38.388' 

12 

2  PUVs 

2 

30°  18.483' 

8 6U  07.493' 

12 

pressure 

3 

29°  59.010' 

85u  34.990' 

17 

pressure 

4 

29°  32.710' 

85°  17.237' 

15 

pressure 

5 

29°  46.247' 

84°  28.861' 

10 

2  PUVs 

pressure 

6 

28°  57.277' 

85°  01.971' 

47 

pressure 

7 

29°  58.053' 

86°  35.724' 

83 

pressure 

8 

30°  13.524' 

8 6U  37.688' 

28 

pressure 

9 

30°  22.908' 

8 6U  38.374' 

18 

wavebuoy 

With  the  exception  of  the  wavebuoy,  all  moorings 
consist  of  a  fiberglass  tripods  (Sea-Spider)  on  which 
instruments  (PUV  or  pressure  sensor)  are  mounted.  Also 
attached  to  the  Sea-Spider  is  an  acoustic  release  with  a 
float  assembly,  which  is  used  during  recovery  of  the 
platform  and  instrument. 
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Figure 
sites . 


2 .  Geographical  location  of  the  nine  deployment 
Contours  indicate  the  Continental  Shelf  bathymetry 
at  20m  intervals. 


Figure  3.  Sea  spider  tripod  with  a  PUV  sensor  in  the 
center.  Mounted  on  the  side  is  an  assembly  with  an  acoustic 

release,  used  for  recovery. 
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Figure  4.  Sea  spider  tripod  with  a  pressure  sensor  in  the 
center  and  dual  acoustic  release  recovering  systems. 

High  quality  data  were  collected  at  six  sites:  three 
to  seven  and  nine.  Instrument  platforms  at  sites  one,  two 
and  eight  were  lost  in  the  extreme  wave  conditions  of 
Hurricane  Ivan. 

C .  SATELLITE  DATA 

The  Navy' s  GEOSAT  Follow-On  (GFO)  satellite  was 
launched  on  February  10,  1998,  following  the  previous 

GEOSAT  that  was  disabled  in  1990. 

GFO  is  a  mini  satellite  designed  and  built  by  Ball 
Aerospace  &  Technology  Corporation  (BATC)  of  Boulder, 
Colorado.  Its  mission  is  managed  by  the 

Meteorological/Oceanographic  (METOC)  Program  Office  of 
Space  and  Naval  Warfare  Systems  Command  (SPAWAR)  in  San 

Diego,  CA,  and  its  major  objective  is  to  provide 
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operational  altimetry  data  for  the  US  Navy  as  well  as  for 
NOAA,  NASA,  and  University  ocean  science  and  ocean 
monitoring . 

A  13.5-GHz  radar  altimeter  is  the  primary  payload, 
providing  wave-height  measurements  to  an  accuracy  of  3.5cm. 
The  altimeter  tracks  the  radar  time  delay,  with  a  closed 
loop  adaptive  algorithm.  If  the  altimeter  transmitted  a 
pulse  and  waited  for  the  return  from  the  surface,  the 
resulting  time  delay  would  be  a  direct  measure  of  the  range 
to  the  surface.  In  order  to  increase  the  ocean  sampling 
rate  and  reduce  height  noise  (bad  measurements) ,  the  radar 
has  five  pulses  in  the  air  at  any  given  time,  and  measures 
the  delay  between  the  most  recent  pulse  generated  and  the 
next  pulse  received  (generated  5  pulse  intervals  ago)  and 
the  true  range  is  computed  by  the  Sensor  Data  Record  (SDR) 
software  on  the  ground  at  the  Payload  Operations  Center 
(POC) ,  co-located  with  the  Altimeter  Data  Fusion  Center 
(ADFC)  at  Naval  Oceanographic  Office  (NAVOCEANO)  at 
Mississippi  (http : //Gf o/Data_val/Cal_f ormats/time_mgt . htm) . 
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III.  WIND  FIELD 


Hurricane  Ivan  was  a  standard  Cape  Verde (1)  hurricane 
that  reached  Category  5  strength  three  times  on  the  Saffir- 
Simpson  Hurricane  Scale (2)  (SSHS)  . 


A.  SYNOPTIC  HISTORY 

Hurricane  Ivan  originated  from  a  large  tropical  wave 
off  the  African  West  Coast  on  August  31,  2004.  Despite  a 

relatively  low  latitude  (9.7°N),  development  continued  and 
it  is  estimated  that  the  cyclone  became  Tropical  Storm  Ivan 
just  12  hours  later  at  0600UTC  September  3  (Stewart, 
2005)  <3)  . 

After  emerging  over  the  southern  Gulf  of  Mexico  early 
on  September  14,  Ivan  turned  north-northwestward  and  then 
northward.  As  Ivan  was  getting  close  to  the  northern  Gulf 
coast,  the  upper-level  wind  flow  ahead  of  the  trough  became 
more  westerly  and  strengthened  to  more  than  30kt,  which 
increased  the  shear  and  advected  dry  air  into  the  inner 
core  region.  The  presence  of  cooler  shelf  water  just 
offshore  weakened  Ivan  before  making  landfall  as  a  105kt 


1-  Cape  Verde  type  hurricanes  are  those  Atlantic  basin  tropical 
cyclones  that  develop  into  tropical  storms  fairly  close  (<1000  km)  to 
the  Cape  Verde  Islands  and  then  become  hurricanes  before  reaching  the 
Caribbean . 

2-  The  Saf f ir-Simpson  Hurricane  Scale  is  a  1-5  rating  based  on  the 
hurricane's  present  intensity.  This  is  used  to  give  an  estimate  of  the 
potential  property  damage  and  flooding  expected  along  the  coast  from  a 
hurricane  landfall.  Wind  speed  is  the  determining  factor  in  the  scale, 
as  storm  surge  values  are  highly  dependent  on  the  slope  of  the 
continental  shelf  in  the  landfall  region. 

3-  For  more  details  concerning  hurricane  Ivan  see  the  following  report: 
Stewart,  Stacy  R.,  Hurricane  Ivan  2-24  September  2004  -  Tropical 
Cyclone  Report,  National  Hurricane  Center,  2005. 

13 


hurricane  (category  3  on  the  SSHS)  at  approximately  0650UTC 
September  16,  just  west  of  Gulf  Shores,  Alabama. 

By  this  time,  the  eye  diameter  had  increased  to  40- 
50NM  (Figure. 5),  which  resulted  in  some  of  the  strongest 
winds  occurring  over  a  narrow  area  near  the  southern 
Alabama-western  Florida  panhandle  border. 

B .  METEOROLOGICAL  MEASUREMENTS 

Observations  during  Ivan  include  satellite-based 
Dvorak14’  technique  intensity  estimates  from  the  Tropical 
Analysis  and  Forecast  Branch  (TAFB) ,  the  Satellite  Analysis 
Branch  (SAB)  and  the  U.  S.  Air  Force  Weather  Agency  (AFWA)  , 
as  well  as  flight-level  and  dropwindsonde  observations  from 
flights  of  the  53rd  Weather  Reconnaissance  Squadron  of  the 
U.  S.  Air  Force  Reserve  Command  (AFRES) . 


Figure  5.  Satellite  picture  of  Hurricane  Ivan  at  16SEP0000Z 

(from  www.noaanews.noaa.gov) 


4-  The  Dvorak  technique  is  a  methodology  to  get  estimates  of  tropical 
cyclone  intensity  from  satellite  pictures.  Vern  Dvorak  developed  the 
scheme  using  a  pattern  recognition  decision  tree  in  the  early  1970s. 
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Microwave  satellite  imagery  from  NOAA  polar-orbiting 
satellites,  the  NASA  Tropical  Rainfall  Measuring  Mission 
(TRMM) ,  the  NASA  QuikSCAT,  and  Defense  Meteorological 
Satellite  Program  (DMSP)  satellites  were  also  helpful  in 
monitoring  Ivan. 

Loaded  with  sophisticated  sensors  and  weather 
measuring  devices,  the  different  aircrafts  were  launched 
every  six  hours  for  missions  that  typically  lasted  eight  to 
twelve  hours.  Aircrews  usually  fly  into  the  building  storm 
at  about  10,000  feet  above  the  ocean's  surface.  During  the 
missions,  the  aircrafts  crisscross  the  hurricane  and 
release  small  "dropsonde"  canisters.  Dropped  by  parachute, 
these  free-floating  sensors  provide  very  accurate 
measurements  of  the  storm's  location  and  intensity.  The 
canisters  relay  barometric  pressure,  wind  speed  and 
direction  and  other  details  to  the  aircraft  during  their 
descent  until  they  hit  water.  After  checking  the  data,  the 
crews  forward  the  information  directly  to  the  National 
Hurricane  Center. 

C.  WINDS  FIELD 

In  order  to  analyze  the  surface  wind  field,  it  is 
fundamental  that  all  the  information  conform  to  a  common 
height.  To  achieve  that,  a  complex  five-step  process  was 
developed  (Powell,  1996)  that  is  briefly  reviewed  here. 
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Step  1  -  Compute  the  mean  (10  min  average) 


The  surface  boundary 
steps  2  to  4  requires  a 
mean  is  used  because  it  i 
mean  winds . 


layer  similarity  theor 
mean  wind  computation, 
s  more  stable  than  the 


y  used  in 
A  10  min 
raw  1  min 


Step  2 


Determination  of 


observation 


exposure 


type 


Figure  6. 


Flow  chart  of  procedures  used 
field  (Powell,  1996) 


to  achieve  a  wind 


In  this  step  it  is  decided  what  type  of  fetch  is  used. 
The  fetch  depends  upon  whether  the  wind  direction  is 
onshore,  offshore  or  alongshore.  The  land  exposure  type  is 
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important  because  wind  blowing  from  sea  to  land  decelerates 
and  forms  an  internal  boundary  layer  (IBL)  whereas  wind 
blowing  from  land  to  sea  accelerates. 

Step  3  -  Adjust  to  the  reference  height  (10  meters) 

The  World  Meteorological  Organization  (WMO)  recommends 
the  wind  measurement  height  to  be  10  meters  above  ground 
level.  The  method  used  to  adjust  wind  speeds  to  that  level 
is  described  by  Liu  et  al .  (1979)  .  The  method  iteratively 
solves  eguations  involving  the  air-sea  exchanges  of 
momentum,  heat,  and  water  vapor  to  arrive  at  the  wind  speed 
profile  in  the  lower  atmospheric  boundary  layer.  Air  and 
water  temperature,  wind  speed  at  the  height  of  observation, 
and  relative  humidity  are  required  inputs  to  the  algorithm. 
Since  relative  humidity  isn't  usually  measured,  it  is 
assumed  to  be  85%. 

Step  4  -  Adjust  to  the  reference  exposure 

In  order  to  simplify  the  analysis,  all  the  input  data 
over  land  is  converted  to  an  equivalent  open  ocean  exposure 
environment  with  a  roughness  length  of  0.03m. 

Step  5  -  Adjust  to  the  maximum  sustained  wind  over  the 
averaging  period 

After  adjustment  to  10m  and  adjustment  to  open  ocean 
for  land  exposures,  observations  with  averaging  times 
greater  than  2  min  are  interpolated  to  1  min  intervals. 


17 


After  completion  of  these  standardization  steps,  the 
data  are  processed  in  an  objective  analysis  scheme,  which 
interpolates  data  from  irregularly  spaced  locations  to  a 
fixed  grid.  A  filtering  procedure  is  applied  in  which 
stations  surrounding  a  nearby  grid  point  receive  a  relative 
weight,  typically  based  on  the  distance  they  lie  from  the 
grid  point.  In  this  way,  the  grid  point  value  is  not 
representative  of  a  single  station,  but  is  instead  a  best 
fit  to  all  the  surrounding  data.  The  effect  of  considering 
several  stations  leads  to  smoothing  of  the  data.  The  final 
grid  point  value  will  range  somewhere  between  the  maximum 
and  minimum  values  of  the  stations  used. 

The  Atlantic  Oceanographic  and  Meteorological 
Laboratory  (AOML)  generated  a  final  gridded  H-Wind  field 
for  hurricane  Ivan  with  spatial  and  temporal  resolutions  of 
0.05°  (~6Km)  and  3  hours  respectively  that  was  used  in  this 
thesis.  The  mobile  grid  is  centered  on  the  eye  of  the 
storm.  Example  wind  fields  at  three  time  steps  are  shown  in 
Figure  7  . 
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WaveWatch  III  Wind  field  (m/s)  09-13-04  00  UTC 


Longitude  (DEG) 
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WaveWatch  III  Wind  field  (m/s)  09-16-04  09  UTC 


-98  -96  -94  -92  -90  -88  -86  -84  -82 

Longitude  (DEG) 


Figure  7.  Example  AOML  H-Wind  field  snapshots  during  Ivan. 
The  wind  field  is  gridded  on  a  moving  "square"  of 
approximately  lOOOKm  by  lOOOKm  centered  on  the  eye  of  the 

storm 
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IV.  WAVE  MODEL 


A.  BRIEF  DESCRIPTION  OF  WAVEWATCH  III  (WW3) 

WAVEWATCH  III  (WW3)  is  a  full-spectral  third- 
generation  wind  wave  model  developed  at  the  Marine  Modeling 
and  Analysis  Branch  (MMAB)  of  the  Environmental  Modeling 
Center  (EMC)  of  the  National  Centers  for  Environmental 
Prediction  (NCEP) .  It  is  based  on  previous  versions  of  the 
same  model  with  improvements  in  coding  structures, 
numerical  methods  and  physical  parameters.  WAVEWATCH  III 
version  2.22  (Tolman,  2003)  was  used  for  the  hindcasts  in 
this  thesis. 

WAVEWATCH  III  solves  the  spectral  wave  action  balance 
equation  for  wavenumber-direction  density  spectra.  The 
implicit  assumption  of  this  equation  is  that  properties  of 
the  medium  (water  depth  and  current)  as  well  as  the  wave 
field  itself  vary  on  time  and  space  scales  that  are  much 
larger  than  the  variation  scales  of  a  single  wave.  A 
further  constraint  is  that  the  parameterizations  of 
physical  processes  included  in  the  model  do  not  address 
conditions  where  the  waves  are  strongly  depth-limited. 
These  two  basic  assumptions  imply  that  the  model  can 
generally  be  applied  on  spatial  scales  (grid  increments) 
larger  than  1  to  10  km,  and  outside  the  surf  zone 
(http : / /polar . ncep . noaa . gov/waves/wavewatch) . 

The  wave  action  balance  is  given  by: 


=  5  =  +  Snl  +  Sds  +  Sbot 


DN 

Dt 
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(1) 


The  left  hand  side  represents  the  total  derivative  (moving 
with  a  wave  component)  of  the  action  density  spectrum  N  and 
S  stands  for  the  net  effect  of  sources  and  sinks.  This  net 
effect  of  sources  in  WW3  is  represented  as  a  superposition 
of  wind  input  (Sin)  ,  non-linear  transfer  (Sni)  and 
dissipation  (Sds)  •  In  order  to  extend  its  application  to 
shallow  coastal  water,  an  additional  bottom  friction  source 
term  (Sbot)  is  included. 

Modifications  and  empirical  fine-tuning  of  the  source 
terms  has  led  to  several  alternative  formulations.  In  this 
study  the  more  recent  formulation,  proposed  by  Tolman  and 
Chalikov  (1996),  was  used.  This  source  term  package 
consists  on  an  input  source,  and  two  dissipation 
components . 

The  input  and  dissipation  terms,  Sin  and  Sds,  are 
treated  as  a  collective  source  term,  since  they  are  inter¬ 
connected  (Tolman,  2002).  The  parameterization  used  is 
based  on  the  fact  that  the  dissipation  processes  for 
frequencies  at  and  below  the  spectral  peak  are  different 
from  those  occurring  at  high  frequencies.  For  this  reason, 
the  WAVEWATCH  III  dissipation  term  contains  two  different 
components,  one  for  low  frequencies,  SdS,i?  dominated  by 
turbulence  and  another  one  for  high  frequencies,  SdS,h/ 
representing  the  losses  by  whitecapping. 

Nonlinear  wave-wave  interactions,  Sni,  were  modeled 
using  a  Discrete  Interaction  Approximation  (DIA)  method 
(Hasselmann  et  al .  ,  1985). 

The  bottom  friction  parameterization  is  based  on  an 
empirical  relationship  from  the  Joint  North  Sea  Wave 
Project  (JONSWAP)  study. 
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B.  GRID  SETTINGS 

For  this  study  two  grids  were  used.  A  large  scale, 
relatively  coarse  grid  with  0.2°  spatial  resolution  covers 
the  entire  Gulf  of  Mexico  and  surrounding  areas  between 
parallels  15°N  and  31.2°N  and  meridians  100°W  and  79.8°W.  A 
fine  resolution  (0.05°)  grid  was  implemented  for  the 
eastern  part  of  the  Gulf  of  Mexico  between  parallels  21° 
and  31°  North  and  meridians  93°  to  80°  West.  High 
resolution  bathymetry  from  the  National  Geophysical  Data 
Center  (NGDC)  was  used  to  accurately  resolve  coastline  and 
bathymetry  effects  in  the  area  of  the  SAX04  experiment.  The 
spectral  calculations  for  both  grids  were  computed  along  24 
directions  spanning  the  full  circle  (every  15°)  and  30 
frequencies  between  0.0412  Hz  and  0.6530  Hz. 

The  wind  field  was  interpolated  in  order  to  match  the 
wave  model  resolution  while  the  bathymetry  was  rearranged. 
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V.  WAVE  MODEL  -  DATA  COMPARISONS 


In  order  to  raise  our  understanding  of  the  variability 
in  wave  conditions  observed  during  Ivan  as  well  as  to  test 
the  accuracy  of  the  WaveWatchlll  model,  several  model  runs 
were  made  and  the  outputs  were  compared  with  the  different 
sets  of  in  situ  data  that  were  collected  by  the  NDBC  buoy 
network  and  in  the  SAX04  experiment,  as  well  as  the  GFO 
satellite  altimeter  data.  Basically,  comparisons  were 
divided  in  two  areas;  Deep  and  Shallow  water.  The  deep 
water  area  covered  the  entire  Gulf  of  Mexico  with  a 
relatively  coarse,  0.2°  degree,  spatial  resolution.  Here 
the  wave-bottom  interactions  are  non-existent,  so  there  is 
no  gain  in  increasing  spatial  resolution.  However,  in  the 
shallow  water  area  of  the  SAX04  Experiment  wave-bottom 
interactions  become  of  main  interest,  so,  a  finer 
resolution  (0.05°  degree)  grid  covering  the  eastern  side  of 
the  Gulf  was  used  to  resolve  the  bathymetry  and  the 
coastline  as  well. 

Dissipation  by  bottom  friction  is  a  very  important 
mechanism  when  comparing  data  in  shallow  water.  In  order  to 
understand  how  much  energy  is  dissipated  by  this  process 
and  how  well  WW3  accounts  for  it,  model  runs  with  and 
without  dissipation  were  compared. 
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Figure  8.  The  left  column  shows  the  synoptic  wind  field  and 
the  right  column  the  corresponding  wave  field  for  Sep.  13 
00:00  UTC,  14  00:00  UTC,  15  00:00  UTC,  16  00:00  UTC  and  16 
09:00  UTC  (landfall) .  The  white  line  indicates  hurricane 
track  with  a  red  filled  circle  at  the  location  of  the 
hurricane  eye.  The  white  arrows  indicate  the  wind  (left 
panels)  and  wave  (right  panels)  directions.  Colors 
correspond  to  wind  speed  (left  panels)  and  significant  wave 

height  (right  panels) . 
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A.  DEEP  WATER 

Before  describing  the  deep  water  model-data 
comparisons,  it  is  useful  to  give  a  synoptic  description  of 
the  predicted  evolution  (Figure  8) .  Since  the  AOML  H-winds 
grid  only  covers  a  square  area  of  approximately  1000  Km  by 
1000  Km,  the  forcing  that  generates  waves  in  WW3  is 
confined  to  this  area. 

The  Gulf  of  Mexico  is  practically  isolated  from  the 
Atlantic  Ocean  and  from  the  Caribbean  Sea  as  well.  The  only 
significant  opening  where  swell  can  enter  the  Gulf  is  the 
Yucatan  Channel  (around  200  Km  wide) .  No  significant  wave 
activity  is  predicted  in  the  Gulf  until  the  storm  covers 
the  Yucatan  Channel  (Figure  8) . 

After  this  point,  waves  have  the  entire  Gulf  to  grow 
and  expand.  In  less  than  24  hours,  waves  reach  a  height  of 
20  meters.  As  the  storm  reaches  the  northern  half  of  the 
Gulf,  a  marked  asymmetry  between  the  wave  field  on  the 
right  and  left  side  of  the  hurricane  track  develops.  On  the 
right  side  (East)  of  the  hurricane  track,  the  waves  are 
propagating  in  the  same  direction  as  the  storm,  allowing 
the  waves  to  spend  more  time  under  the  storm  influence 
resulting  in  higher  waves.  On  the  other  hand,  on  the  left 
side,  waves  are  propagating  away  from  the  storm  resulting 
in  smaller  waves. 

As  described  in  chapter  2,  seven  NDBC  buoys  were 
chosen  for  comparisons.  Model-data  comparisons  of  the  bulk 
parameters.  Significant  Wave  Height  (SWH) ,  Mean  Wave 
Direction  (MWD)  and  Peak  Period  (Tp)  are  shown  in  Figures 
9-11 . 
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Figure  9.  Comparison  of  observed  (blue  curve)  and  predicted 
(red  curve)  significant  wave  heights  at  the  seven  NDBC 

buoys . 

Since  the  seven  buoys  are  all  located  in  the  northern 
half  of  the  Gulf,  no  significant  wave  activity  (SWH  close 
to  zero)  is  predicted  on  the  first  day  (September  13)  when 
the  storm  enters  the  Gulf  (Figure  9)  .  The  observed  SWH 
reflect  local  atmospheric  forcing  that  is  not  accounted  for 
in  the  WW3  predictions.  There  is  a  notable  time  difference 
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between  when  maximum  wave  heights  were  observed  versus  the 
WW3  prediction.  After  the  storm  passes  by  and  the  wave 
heights  start  to  decrease,  WW3  predictions  agree  reasonably 
well  with  the  observations.  Buoy  42040  was  damaged  during 
the  passage  of  the  hurricane  and  did  not  report  data  during 
the  waning  stage  of  the  storm. 

The  predicted  mean  wave  directions  are  in  reasonable 
agreement  with  the  observations  (Figure  10) .  Predictions  at 
buoys  42001  and  42003  capture  the  observed  direction 
inversions  as  the  storm  passes  by,  with  waves  first 
arriving  from  East  quadrants  and  then  turning  to  North  or 
Northwest . 

To  describe  the  directional  wave  field  in  more  detail, 
predictions  of  the  frequency-directional  wave  spectra 
during  and  after  the  storm  passage  are  compared  in  Figures 
11  and  12.  At  station  42001,  on  the  West  side  of  the  storm, 
the  energy  at  lower  frequencies  arrives  from  the  East 
spreading  to  the  northern  quadrants  at  higher  frequencies. 
This  pattern  is  consistent  with  the  cyclonic  wind 
circulation  (Figure  11-a) .  Whereas  the  lower  frequency 
swell  responds  to  forcing  over  several  days  from  the 
predominant  easterly  winds,  the  high  frequency  part  of  the 
spectrum  quickly  adjusts  to  the  local  wind  direction,  which 
is  from  the  North.  Later  on  September  16  (Figure  11-b) ,  the 
storm  had  already  made  landfall,  and  the  wind  has  been 
blowing  from  North  with  a  restricted  fetch.  The  energy  is 
concentrated  at  higher  frequencies  and  arrives  from  the 
northern  quadrants,  with  a  smaller  residual  swell  from  the 
east . 
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Figure  10.  Comparison  of  the  observed  (blue  crosses)  and 
predicted  (red  crosses)  mean  wave  direction  at  six  NDBC 

buoys . 

On  the  other  hand,  at  station  42003  located  on  the 
East  side  of  the  storm,  (Figure  12-a)  ,  the  dominant  low 
frequency  wave  energy  during  the  storm  passage  arrives  from 
the  south.  Later,  as  the  storm  made  landfall,  (Figure  12-b) 
the  energy  distribution  is  similar  to  that  predicted  at 
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Frequency  (Hz)  Frequency  (Hz) 


station  42001,  but  with  low  frequency  residual  swell 
arriving  from  the  northwest,  again  consistent  with  the 
cyclonic  circulation. 
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a)  b) 

Figure  11.  Two-dimensional  frequency-directional  spectra 
predicted  by  WW3  at  station  42001.  a)  15SEP  09:00UTC;  b) 
16SEP  15:00  UTC.  Contours  indicate  energy  levels  on  a  log 
scale  at  lO17-  intervals  relative  to  the  maximum  energy. 


Figure  12.  Two-dimensional  frequency-directional  spectra 
predicted  by  WW3  at  station  42003.  a)  15  SEP  09:00UTC;  b) 
16SEP  15:00  UTC.  Contours  as  in  previous  figure. 
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igure  13.  Comparison  of  observed  (blue  curve)  and  predicted 
(red  curve)  peak  period  at  the  seven  NDBC  buoys. 


Regarding  the  peak  period  (Figure  11),  on  September  13 

and  14  there  are  large  differences  between  the  model 

predictions  and  observations.  The  model  predicts  very  low 

frequency  (and  low  energy)  waves  that  are  the  forerunners 

of  Ivan  whereas  the  observed  wave  field  at  the  NDBC  buoys 

is  still  dominated  by  higher-frequency  locally  generated 

seas  (not  accounted  for  in  the  WW3  hindcast) .  The  agreement 
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between  observations  and  predictions  improves  dramatically 
on  the  latter  two  days  when  the  large  waves  from  Ivan  reach 
the  buoys . 

The  satellite  altimeter  measurements  allow  for 
additional  WW3  model  verification  in  regions  of  the  Gulf  of 
Mexico  where  there  are  no  NDBC  buoys.  Observed  and 
predicted  significant  wave  height  along  seven  satellite 
tracks  that  overlap  the  model  domain  (Figure  14)  are 
compared  in  Figure  15.  Results  confirm  the  earlier  noted 
model  tendency  to  underpredict  SWH. 


Figure  14.  GEOSAT  FOLLOW  ON  (GFO)  satellite  tracks  - 

September  13-16,  2004 
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Figure  15.  Comparisons  between  GEOSAT  observations  and 
WaveWatch  III  predictions  along  the  tracks  shown  in 

Figure  14. 

Red  line  -  WW3  predictions 
Blue  line  -  GEOSAT  measurements 


B.  FLORIDA  CONTINENTAL  SHELF 

To  resolve  coastal  boundary  and  bathymetry  effects,  a 
high  resolution  (0.05°  degree)  grid  was  implemented  for  the 
Eastern  part  of  the  Gulf  including  the  Florida  panhandle 
shelf  where  the  SAX04  measurements  were  colleted.  To 
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investigate  the  importance  of  dissipation  (bottom  friction) 
on  the  continental  shelf  and  evaluate  the  empirical  JONSWAP 
parameterization  used  in  WW3,  model  runs  were  performed 
both  with  and  without  bottom  friction. 

The  area  where  Hurricane  Ivan  made  landfall  features  a 
complex  coastline.  From  Venice,  LA  to  Panama  City,  FL,  the 
coastline  is  highly  irregular  with  a  chain  of  barrier 
islands  and  a  complex  inner  shelf.  The  barrier  islands 
south  of  Pensacola  and  Fort  Walton  Beach  form  a  perfect 
natural  wave  barrier.  Since  nearshore  processes  (i.e.  storm 
surge  and  surf  zone  wave  breaking)  are  not  accounted  for  in 
WW3,  the  model  computations  were  restricted  to  areas  deeper 
than  5  meters . 


WaveWatch  III  Wave  field  (m)  09-15-04  21  UTC 
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WaveWatch  III  Wave  field  (m)  09-1 6-04  03  UTC 
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WaveWatch  III  Wave  field  (m)  09-1 6-04  06  UTC 


Figure  16.  WW3  Wave  field  (SWH)  predictions  for  September  15 
21:00UTC  and  September  16,  03:00  and  06:00  UTC 
The  Hurricane  track  is  represented  by  the  white  line 
and  the  red  spot  indicates  the  eye  of  the  storm. 


The  Highest  SWH  is  predicted  to  the  east  and  ahead  of 
the  storm  center.  As  discussed  earlier,  on  the  east  side, 
waves  propagate  along  with  the  storm  thus  increasing  the 
"dwell  time"  when  waves  are  under  the  influence  of 
hurricane  force  winds.  Since  group  speed  of  the  dominant 
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waves  is  large  than  the  speed  with  which  the  storm 
advances,  the  largest  waves  are  found  slightly  ahead  of  the 
storm  center.  On  the  west  side,  waves  propagate  away  from 
the  storm,  minimizing  the  time  that  they  are  under  the 
storm  influence,  resulting  in  lower  wave  heights. 

Comparing  these  high-resolution  synoptic  fields  with 
the  earlier  discussed  coarser  resolution  run,  the 
differences  between  the  predicted  mean  wave  directions  are 
surprisingly  small  even  in  the  shallow  shelf  areas.  These 
results  indicating  that  refraction  effects  are  weak.  On  the 
right  side  of  the  hurricane  track,  the  unlimited  fetch 
drives  long  period  waves  that  feel  the  bottom  across  the 
entire  shelf,  but  these  waves  propagate  northward  almost 
perpendicular  to  the  bathymetry,  minimizing  the  refraction 
effects.  On  the  other  hand  on  the  left  side  of  the  storm, 
waves  propagate  southeastward,  and  since  the  wind  is 
blowing  offshore  the  fetch  is  very  limited,  resulting  in 
high  frequency  waves  that  also  are  not  strongly  affected  by 
refraction  (Figure  17) . 


WaveWatch  III  Peak  Period  field  (s)  09-16-04  06  UTC 


Figure  17.  WW3  Peak  Period  field.  The  Hurricane  track  is 
represented  by  the  white  line  while  the  red  spot  indicates 

the  eye  of  the  storm. 
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Comparisons  WW3  -  Pressure  Sensor  site  7 


Time  (decimal  days) 


Comparisons  WW3  -  Pressure  Sensor  site  9 
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Figure  18.  Comparison  between  WW3  predictions  and 
observations  at  SAX04  sites  7  and  9.  Blue  line  -  WW3 
predictions  with  dissipation 
Red  line  -  SAX04  instruments  observations 
Green  line  -  WW3  predictions  without  dissipation 
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Figure  19.  Comparison  between  WW3  predictions  and 
observations  at  SAX04  site  5. 

Blue  line  -  WW3  predictions  with  dissipation 
Red  line  -  SAX04  instrument  observations 
Green  line  -  WW3  predictions  without  dissipation 

At  SAX04  site  7  (83  meters  depth  -  Figure  18)  the  WW3 

predictions  with  and  without  dissipation  are  practically 
identical,  indicating  that  bottom  friction  in  these 
relatively  deep  regions  is  negligible.  On  the  other  hand, 
at  site  9,  WW3  predictions  without  dissipation  match  the 
actual  in  situ  observations  almost  perfectly  whereas 
predictions  with  dissipation  are  biased  low  by  about  40%. 
Although  it  is  possible  that  WW3  is  dissipating  too  much 
energy,  the  agreement  without  dissipation  must  be 
fortuitous,  since  the  large  waves  propagating  through  this 
shallow  region  are  approaching  depth-limited  wave  breaking. 
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Overall,  these  comparisons  indicate  that  the  wave  decay  by 
bottom  friction  is  relatively  weak  across  this  narrow 
section  of  the  continental  shelf. 

At  site  5,  inshore  of  the  wide  shelf  of  Apalachee  Bay, 
WW3  predicts  large  differences  in  wave  heights  with  and 
without  dissipation,  consistent  with  the  strong  dissipation 
expected  across  a  wide  shelf.  The  model  without  dissipation 
overpredicts  the  SWH  by  almost  a  factor  of  2  whereas  the 
model  with  dissipation  underpredicts  the  SWH  by  about  30%, 
suggesting  that  the  default  setting  for  dissipation  induced 
by  bottom  friction  in  WW3  is  too  high  for  these  conditions. 
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VI.  CONCLUSIONS 


The  objective  of  this  study  was  to  test  the 
performance  of  the  operational  wave  prediction  model 
WaveWatch  III  in  regional  applications  under  extreme 
meteorological  conditions,  such  as  hurricanes.  Model 
hindcasts  were  performed  for  category  V  Hurricane  Ivan 
which  made  landfall  on  the  Alabama  coast  in  2004.  Model 
predictions  are  compared  with  measurements  collected  in  the 
ONR  Sediment  Acoustics  Experiment  SAX04  and  data  from  NDBC 
buoys  in  the  Gulf  of  Mexico  and  the  Geosat  Follow  On 
Satellite . 

Three  model  runs  were  made.  One  run,  covering  the 
entire  Gulf  of  Mexico,  with  0.2°  degree  spatial  resolution 
and  two  runs,  for  the  Eastern  part  of  the  Gulf,  with  0.05° 
degree  resolution,  one  with  and  one  without  dissipation. 
The  model  was  forced  with  AOML  H-winds . 

Comparisons  with  seven  NDBC  buoys,  13  SAX04 
instruments  and  GFO  altimeter  data  were  made  for  the  period 
of  September  13  to  September  16,  2004,  while  Hurricane  Ivan 
crossed  through  the  Gulf  of  Mexico. 

In  deeper  areas,  analysis  reveals  that  WaveWatchlll 
has  a  tendency  to  underpredict  wave  heights  and  to  be 
consistently  late  (about  four  to  six  hours)  in  predicting 
the  arrival  of  the  dominant  waves.  On  the  other  hand,  the 
predicted  decaying  phase  of  the  wave  evolution  matches  the 
observations  fairly  well.  In  shallow  areas,  the  model 
predicts  too  much  dissipation. 

The  model  predicts  large  differences  between  wave 
conditions  on  the  left  and  right  sides  of  the  hurricane 
path  owing  to  the  difference  in  "dwell  time"  between  waves 
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propagating  against  and  with  the  storm.  The  reduced  fetch 
on  the  left  side  of  the  storm  limits  the  wave  growth  during 
landfall . 

Overall,  the  WaveWatchlll  hindcasts  are  in  reasonable 
agreement  with  the  observations. 
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